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Precise measurements of the baryon acoustic oscillation (BAO) scale as a standard ruler in the 
clustering pattern of large-scale structure is a central goal of current and future galaxy surveys. 

The BAO peak may be sharpened using the technique of density-field reconstruction, in which the 
bulk displacements of galaxies are estimated using a Zel’dovitch approximation. We use numeri¬ 
cal simulations to demonstrate how the accuracy of this approximation depends strongly on local 
environment, and how this information may be used to construct an improved BAO measurement 
through environmental re-weighting and using higher-order perturbation theory. We outline further 
applications of the displacement field for testing cosmological models. 


I. INTRODUCTION 

The large-scale structure of the Universe provides one 
of our most powerful tests of the cosmological model, 
encoding a wealth of information about the expansion 
history of the Universe (imprinted as a standard ruler in 
baryon acoustic oscillations) and its gravitational physics 
(inferred from the growth of structure with time). A cen¬ 
tral goal of modern cosmology is to map out this struc¬ 
ture through galaxy redshift surveys, and to develop ac¬ 
curate models to link these observations to theory. 

In the modelling of cosmological fluid dynamics, a par¬ 
ticularly important role is played by the vector displace¬ 
ment field (q, t ), which specifies the trajectory of fluid 
elements through space (here denoted by a Lagrangian 
co-ordinate q) and time t. In modelling approaches such 
as Lagrangian Perturbation Theory, perturbative solu¬ 
tions may be formulated for ’S' and used to construct 
models for the statistics of the observed density field 

(e.g.EHi) 

The displacement field has assumed particular impor¬ 
tance in “density-field reconstruction” (@), a technique 
used to sharpen measurements of the baryon acoustic 
peak in the galaxy correlation function. Measurements 
of the baryon acoustic peak, a preferred clustering scale 
imprinted in the distribution of photons and baryons by 
the propagation of sound waves in the relativistic plasma 
of the early Universe (d, 0), have assumed particular 
importance in recent years as a robust standard ruler 
to map out the cosmic expansion history ([1, 0). How¬ 
ever, the pristine sound-horizon scale imprinted in the 
high-redshift matter distribution is “blurred” by the bulk 
displacement of galaxies from their initial positions f[lB- 
13). Reconstruction seeks to use the observed density 
field to estimate these displacements, for example using 
the Zel’dovitch approximation ([la,[l3|) (ZA), and hence 
retract galaxies to their near-original positions in the den¬ 
sity field, restoring the sharp preferred separation. The 
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technique has been successfully applied to datasets from 
the Sloan Digital Sky Surve y the Baryon Oscilla¬ 

tion Spectroscopic Survey ([l3) and the WiggleZ Dark 
Energy Survey ( 0 )- 

A series of studies has considered the theory and 
implementation of density-field reconstruction. It may 
be formulated within Lagrangian Perturbation Theory 
(@) mi), to which the ZA is the lowest-order contribu¬ 
tion, or by alternati ye p erturbative schemes ([12, Ull), 
or in Fourier space ([l^)- The method is robust to the 
treatment of the bias of the tracer or redshift-space ef¬ 
fects (Pmiii), although room for improvement in the 
algorithm certainly exists (in, Hi). 

In this paper we undertake a forensic study of the per¬ 
formance of the current reconstruction algorithm through 
comparison with the exact displacement field derived 
from particle-tracking in N-body simulations. First, we 
study the accuracy of the ZA using the initial matter 
velocity field with different hypotheses for the smooth¬ 
ing procedure and its sensitivity to the local density. 
This fiducial investigation is not affected by complica¬ 
tions due to galaxy bias and non-linear effects on the 
estimated density field. Using this to create an unbiased 
understanding of how proto-halo displacement is opti¬ 
mally modelled, we extend our results to the reconstruc¬ 
tion of the displacement from the non-linear low-redshift 
halo density field. In particular, we study the error in 
the estimated displacements as a function of galaxy halo 
mass, redshift, filter and smoothing scale for estimating 
the density field. We also also study how this error de¬ 
pends on the order of perturbation theory applied, and 
local environment. We discuss the consequences of these 
results for the optimal measurement of the distance scale 
and the use of the displacement field to test cosmological 
models. For the purposes of this investigation, we neglect 
redshift-space distortions in the galaxy co-ordinates. 

Our paper is structured as follows: in Section II we re¬ 
view the connection between the displacement and den¬ 
sity fields within Lagrangian Perturbation Theory. We 
describe our methods for recovering the exact displace¬ 
ment field from N-body simulations, and the estimated 
displacement field from the initial simulation velocity 
field. In Section III we examine the performance of the 
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density-field reconstruction algorithm, in particular fo¬ 
cussing on the degradation of this performance with in¬ 
creasing local density. In Section IV we suggest how the 
measured post-reconstruction correlation function can be 
“re-weighted” by environment to maximize the sharpness 
of the baryon acoustic peak, and we consider the level of 
improvement which may be obtained in the fitted dis¬ 
tance scale. In Section V we discuss the implications of 
our findings for the use of density-field reconstruction, 
and the displacement field itself, for testing cosmological 
models. 


II. DISPLACEMENT FIELD: THEORY AND 
MEASUREMENT 

A. The Zel’dovitch and 2LPT approximations 

In this section we briefly review Lagrangian Pertur¬ 
bation Theory (LPT) at first (i.e., ZA) and second or¬ 
der (i.e., 2LPT). The displacement vector field links 
an initial (Lagrangian) position q of a mass element M, 
smoothed on a scale Rs{M), to the current (Eulerian) 
position, 


x(i?s) = q(i?s) +'®'(q, Als). (1) 

In the ZA, fit is determined by the gradient of the grav¬ 
itational potential V$, which can be expressed in terms 
of the linearly extrapolated density field. Considering 
only the scalar contribution at first and second order 
in the matter density field, we have (i) 

^ ~ ^(1) _|_ iy[r(2) 

( 2 ) 

where D is the linear growth factor and the gravitational 
potential can be solved using the following Poisson equa¬ 
tions: 


v2</>(i)(i?s) = 5o(i?s) 
v2</,(2)(i?s) = -i52(<(>('^) 

where Sq{Rs) is the linear density field smoothed on scale 
Rs and 


a2(<^(i)) = ^[(V.,</.«)2-(V2</>(i))2 


(4) 


The first order displacement (^^^ term) is the well-known 
ZA which can be expressed as a function of the initial 
peculiar velocity 

v{Rs) = adtD{z)dD^^^\ (5) 


where a is the cosmic scale factor. The second order term 
has the label . The variance of the displacement field 


in the ZA is directly proportional to the initial matter 
power spectrum Q such that (e.g. for the x-component) 

(^'W(i?s)2) W\k,Rs)PL{k,Zi)dk, (6) 

where Pl is the linear matter power spectrum at initial 
redshift Zi and W is the Fourier transform of the smooth¬ 
ing function. The higher orders of the displacement field 
can be expressed as higher-dimensional integrals [^. 

In order to determine the accuracy of the ZA in various 
scenarios, we test two cases: 

• Case I: we use the proto-halo velocities to predict 
halo positions at z = 0 assuming the ZA, and com¬ 
pare to the exact simulation displacement field (see 
Sec II B,C). 

• Case 2: we use the halo density field at 2 = 0, and 
second order LPT, to reconstruct the displacement 
— fP of halos, and compare to the exact simulation 
displacement field (see Sec III). 

The first case extends the work of Q (for the un¬ 
smoothed dark matter density field) to the displacement 
of proto-halos. In particular, this analysis provides a 
physical insight of how proto-halos are displaced from 
their positions in the initial matter density field through 
cosmic time up to z = 0. We also extend these studies to 
investigate the sensitivity of the result to the smoothing 
scale and local environment, yielding useful comparisons 
with case 2. 

The second case has a powerful application in the re¬ 
construction of the baryon acoustic peak. Previous stud¬ 
ies have applied this method using first order LPT and 
found a low sensitivity to the smoothing scale Rs which 
enters into Eq. (12131) ([H [13). We extend these tests 
to consider the effect of local environment and evaluate 
the gain of adding the 2LPT correction to reconstruct 
the displacement of halos. In particular, we will show 
how the second order correction depends on the choice of 
smoothing scale and environment. 


B. Exact displacement field from N-body 
simulations 

In order to test the accuracy of the displacement field 
in cases 1 and 2, we use the DEUS simulations. These 
simulations were run for several scientific purposes de¬ 
scribed in [33 113, m, ilgj - idlll . The box side of the sim¬ 
ulations is 648Mpc and they contain 1024^ parti¬ 
cles of mass ~ 1.75 x The simulations were 

carried out using the RAMSES code [13 for a ACDM 
model calibrated to the WMAP 5-year cosmological pa¬ 
rameters [ 13 . The halos were identified using the Friend- 
of-Friends (FoF) algorithm with linking length 6 = 0.2. 

In what follows we will use the Lagrangian size of ha¬ 
los Rl where their mass is given by M{Rs) = V{Rl)p, 
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where p is the mean matter density. We will approximate 
the shape of proto-halos to be spherical, with volume 
v{Rl) = 

We measured the position of the center of mass x/ of 
all halos at z = 0. In order to measure the position of 
the proto-halos in the initial conditions, we labelled all 
particles which belong to a halo at z = 0 and compute 
the center of mass Xi for these particles in the initial 
conditions of the simulation. Therefore we measured for 
each halo the exact displacement = Xf — Xi. 


C. Displacement prediction from initial velocity 
field 

The final position of a proto-halo at z = 0 may be pre¬ 
dicted from the initial velocity field according to Eq.([5|). 
Hence, the displacement of a mass element contained in a 
smoothing radius Rs can be expressed as function of the 
average peculiar velocity initially contained in a patch of 
size Rs, Vi{q,Rs) as: 


'®'«(Q,Z,i?s) 


Vi{Rs)D{z) 

aiH{ai)f{ai)D{zi) 


(7) 


where H{ai) is the Hubble parameter evaluated at the 
initial scale factor and / = din D/d In a is the linear 
growth rate. In this approximation we can use Eq. © to 
displace the proto-halos by measuring the average veloc¬ 
ity Vi around the center of mass q = x^ and within a 
smoothing radius Rg- We investigated the accuracy of 
the ZA in predicting the magnitude of = | ’S' |. 

The smoothing scale Rs is a free parameter in our 
investigation, used in the estimation of the local over¬ 
density field. For Rg —?► oo the average velocity within 
the smoothed region tends to zero, and hence the proto¬ 
halos are not displaced. On the contrary, as Rs —0 we 
obtain the displacement of individual particles, without 
considering that they belong to a halo. 

In Fig.(IT]), we show the measured PDF of S' (black 
squares) for two different masses of halos: the up¬ 
per panel corresponds to halos with Lagrangian radii 
Rl = 5h~^ Mpc (M = while the lower 

panel corresponds to halos with Lagrangian radius Rl = 
8.7 h~^ Mpc (M = As we can see, the 

mean displacement is not strongly sensitive to the mass 
of the proto-halos (varying by only ^ 15%). However, the 
variance of the absolute displacement is reduced for large 
halos. Physically, this is because larger proto-halos are 
less sensitive to the external shear field which contributes 
to the overall dynamics. In Fig.(IT]), we can also see the 
prediction of Eq.© (plotted as histograms) using differ¬ 
ent smoothing scales Rs (see legend) when measuring the 
average peculiar velocities. For Rg < 2 Rl, the absolute 
displacement predicted by Eq. © is approximately inde¬ 
pendent of Rs- For Rs = 4 ,Rl, we begin to smooth 
out relevant density fluctuations and hence underesti¬ 
mate the absolute displacement. In the limit Rg » Rl, 
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FIG. 1: The distribution of Vl' measured in N-body simula¬ 
tions for two different halo masses (black squares). The upper 
panel shows halos with Lagrangian radius Rl = 5h~^ Mpc 
while the bottom panel shows halos with Lagrangian radius 
Rl = 8.7 h~^ Mpc. For each panel, the histograms show the 
Zel’dovitch prediction Eq. © using different smoothing scales 
Rs (see legend). The solid lines correspond to the Maxwell- 
Boltzmann distribution with variance (7%{Rs) computed ac¬ 
cording to Eq. ®. 


the analytical pdf of iF is given by a Maxwell-Boltzmann 
distribution 1|4^L 


n(^',cr^) 




exp 




( 8 ) 
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FIG. 2: Standard deviation of the difference between the x- 
component of the exact displacement and the ZA Eq. 0 , as a 
function of the ratio of the smoothing length and Lagrangian 
size of the halos, Rs/Rl- The red, black and blue crosses 
correspond to halo masses 10^^ ®, 10^®'° and Mq. 


where the variance is 

al,{R) = -^ J W‘^{R,k)PL{k,z)dk. (9) 

The Fourier transform of the filter function W{R,k) is 
chosen to be the Fourier transform of a top-hat filter in 
real space. In Fig. © we can see this analytical prediction 
for R = ARl and R = 2Rl (solid lines). For R = ARl the 
analytic prediction of Eq.® matches the Zel’dovitch pre¬ 
diction of Eq. © with Rs = ARl (comparing the green 
solid line and green histogram). For smaller smoothing 
scales {R = 2Rl), the regime R» Rl is no longer sat¬ 
isfied and we start observing deviations between Eq.® 
and Eq.0 (comparing the blue solid line and blue his¬ 
togram). Interestingly, we observe that Eq.Q qualita¬ 
tively follows the trend for the absolute displacement of 
proto-halos (black squares) once we use R = 2Rl. This 
observation could be used to model the clustering of halos 
with high accuracy, although this goes beyond the scope 
of this paper. 

Second, we consider the dispersion between the exact 
and estimated displacements. Eor each component i, the 
PDE of d'i — follows a Gaussian distribution with 
a mean value equal to zero. In FigI2] we can see the 
standard deviation of the difference between the exact 
and the Zel’dovitch prediction (Eq.0) of 4'a,, for halo 
masses 10^^ ®, 10^^ ° and Mq, as a function of 

the smoothing length Rs (similar results are obtained for 
Afy and 4'^). Independently of the proto-halo masses, the 
optimal smoothing scale to predict the final position of 


FIG. 3: Standard deviation of the difference between the ex¬ 
act displacement and that predicted by Eq. (1^ with Rs = 
Rl in different environments. Colors represent different envi¬ 
ronments. The black square corresponds to the average den¬ 
sity around halos 5{Rs = Rl)- The blue, red and orange 
squares correspond to environments with local overdensity 
5 = {5 — (Ta, J + (Ta, (5 -I- 2as} respectively, where <5 is the mean 
overdensity and us is the standard deviation of the overden¬ 
sity distribution amongst the halos. 


halos is Rs = Rl- At this scale 41 can be determined 
with an accuracy of ^ 0.5 h~^ Mpc, a few per cent of the 
mean displacement. This is an impressive result consid¬ 
ering that we only use the ZA. It implies that if we can 
identify proto-halos in the initial matter density field, we 
can predict their clustering at z = 0 with high accuracy 
in a very short time. 

Finally, we test how sensitive the ZA (Eq.0) is to 
the local environment. Each environment is defined by 
the number density of halos inside a sphere of radius Rs- 
We divide the halos into environment bins by considering 
the probability density function for the number of halos 
at z = 0 within Rs = Rl, in particular the mean 6 
and standard deviation as of this distribution. We then 
define 4 different environments corresponding to halos 
with surrounding densities 5 = {5 — as,5.,5 + as,5 + 2as}- 

In Eig.0 we can see the standard deviation be¬ 
tween the true value and the ZA prediction of Atx (with 
Rs = Rl), as function of the surrounding number 
density of halos {5{Rs)), considering halos with mass 
M = 10^^ h~^MQ, which is equivalent to Rl ~ 3h~^ 
Mpc. The black square corresponds to the mean over¬ 
density environment. As we can see, the accuracy of the 
ZA prediction is sensitive to the environment. Eor over- 
dense environments, the ZA prediction becomes less ac¬ 
curate compared to underdense environments. This sug- 
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FIG. 4: Standard deviation of the difference between the x- 
component of the exact displacement, and the ZA computed 
using the 2 = 0 density field using Eq. (1^ . as a function of 
the smoothing length Rs and the type of filter used. 


gests that the displacement field of proto-halos contains 
non-linear information that can not be encapsulated in 
the simple ZA. An extension of this work would be to 
investigate if the sensitivity to the smoothing scale and 
environment depends on cosmology. 

In the next section, we will repeat this exercise, com¬ 
paring the displacements reconstructed using the halo 
distribution at 2 = 0 to evaluate the density contrast in 
Eq.®, to the exact displacements. 

III. PERFORMANCE OF DENSITY-FIELD 
RECONSTRUCTION 

A. Displacement prediction from density-field 
reconstruction 

We first describe how we predict the displacement field 
from the z = 0 halo distribution using the technique of 
density-field reconstruction. Dividing the volume into a 
3D grid, we first estimated the halo overdensity distribu¬ 
tion at each grid point, using a particular filter (Gaus¬ 
sian or top-hat) and smoothing scale. We then con¬ 
vert the halo overdensity Sh{Rs) to a matter overden¬ 
sity Sm{Rs) = Sh{Rs)/b assuming a linear bias factor 
b, which we determine by comparing the amplitude of 
the large-scale halo power spectrum to the linear-theory 
power spectrum used to generate the simulation. This 
allows us to solve Eq.d^]) using Fast Fourier Transform 
methods, and hence determine the displacement field 
—via Eq.®. 


FIG. 5: Standard deviation of the difference between the x- 
component of the exact displacement, and that computed us¬ 
ing the 2 = 0 density field using Eq. (12131) . as a function of the 
smoothing length Rs using a Gaussian filter. Black squares 
show the first-order (ZA) prediction while black crosses show 
the prediction including the 2LPT correction. 


B. Dependence on smoothing scale and filter 

Considering only the first order of ll' in Ea. (j2l3|) . we 
used several smoothing scales in the range 0 < i?s < 
25 h~^ Mpc for halos with a mass M > 
typical of Luminous Red Galaxies probed by large-scale 
structure surveys such as the Baryon Oscillation Spectro¬ 
scopic Survey ([H). We used both a top-hat filter and 
a Gaussian filter, and we measured the error made in 
the Zel’dovitch prediction compared to the exact values 
for each component as a standard deviation cr('0p^ —'*Ai)- 
The result for the a;-axis is shown in Fig. other axes 
produced similar results. We find that there is again 
an optimal smoothing scale which depends on the filter 
function. The optimal performance of both filters leads 
to a residual error cr ~ 2.5 h~^ Mpc, around 5 times 
poorer than the determination based on the initial proto¬ 
velocities. 

For small smoothing scales the performance of the 
reconstruction method becomes worse, since the num¬ 
ber of neighbours becomes too small to reconstruct the 
density field with sufficient accuracy. Furthermore, the 
higher order correction in Eq.Q becomes non-negligible 
when —>• 0. For large smoothing scales, the degrada¬ 
tion from the optimal performance is a slow function of 
i?s, consistent with the displacement being generated by 
large-wavelength modes that may be successfully recov¬ 
ered even in the presence of significant smoothing (i). 
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FIG. 6: Standard deviation of the difference between the ex¬ 
act displacement and that predicted by Eq. (1^ based on the 
z = 0 halo density field. The squares show the first order 
implementation of Eq. (I2]3t while the crosses implement the 
second order correction. Colors represent different environ¬ 
ments defined in the text (from bottom to top: El, E2, E3, 
E4 and E5). For each color we choose 3 smoothing scales 
(from left to right: Rs = 10, 7.5, 5Mpc. 


These effects are illustrated in Fig.®. Using a Gaus¬ 
sian filter we show the error made in the Ist-order (black 
squares) and 2LPT (black crosses) predictions for differ¬ 
ent smoothing scales. The 2nd-order correction leads to 
a ~ 2% improvement for Rs ~ 8h~^ Mpc and a ~ 12% 
improvement for Rs 5h~^ Mpc. In the next section 
we will emphasize the connection between this result and 
the local environment. 

Overall, the choice of filter (Gaussian vs. top-hat) does 
not affect the performance of the displacement recon¬ 
struction at first-order LPT, once we choose the optimal 
smoothing scale (see Fig.Q). Therefore in what follows 
we will restrict our analysis to a Gaussian filter. 


C. Dependence on local environment 

In this Section we investigate how the accuracy of 
the reconstructed displacement field depends on the lo¬ 
cal environment of each halo. Given a smoothing scale 
Rs, we can compute the number of neighbours for each 
halo and therefore estimate the local density contrast. 
We used the mean local density around halos S, and 
the standard deviation as, to define five bins in envi¬ 
ronment El, E2, E3, E4 and E5 with central values 
(5 = { 5 — cr^, j -I- as, 6 + 2as, 6 + Scr^}. Eor each local 
environment we measured the error in the reconstructed 


displacement at different smoothing scales. 

In Fig.® we show the standard deviation of the differ¬ 
ence between the reconstructed and exact ikxj as function 
of the surrounding smoothed halo density S{Rs)- The 
blue, black, red, orange and green data points correspond 
respectively to environments El, E2, E3, E4, and E5. For 
each environment we plot 3 different squares and 3 differ¬ 
ent crosses which correspond to three different smoothing 
scales Rs = 10,7.5,5/i“^ Mpc, from left to right. The 
squares correspond to the prediction of first-order LPT, 
while the crosses implement the second-order correction. 

In agreement with Fig.®, the accuracy of the predic¬ 
tion is highly sensitive to the environment. Furthermore, 
in Fig.® we can see the effect of the second-order cor¬ 
rection. Unsurprisingly, if we choose a large smoothing 
scale (e.g. 10 h~^ Mpc), the second order correction in 
Eq.® is negligible. Eor a smaller smoothing scale (e.g. 
Rs = 5 h~^ Mpc), the second-order correction becomes 
more important. In fact, the second-order correction con¬ 
tains a term proportional to S^{Rs) and for small values 
of Rs the average local density around halos increases 
{S{Rs = 10) ~ 0.8, S{Rs = 5) ^ 2.8). However we 
can observe that even for Rs = 5 h~^ Mpc, the improve¬ 
ment in the residual from the second-order correction is 
less significant than the difference in residual between 
El and E2. This suggests that the displacement con¬ 
tains non-linear information that can not be described 
by deterministic corrections. 

Einally, for dense environments (e.g 6 > S), a smaller 
smoothing scale leads to a better reconstruction of the 
final halo position predicted by the 2LPT approxima¬ 
tion (Eq® . For halos in a region with the mean density 
or below, reducing the smoothing scale does not partic¬ 
ularly improve the accuracy of the displacement field. 
This convergence of the optimal smoothing scale toward 
Rs —>■ Rl shows the limits of the reconstruction method. 
Reconstruction applied at z = 0 on the non-linear halo 
density field is not equivalent to predicting the final po¬ 
sition of halos from the initial conditions. 


D. Dependence on redshift and halo mass 

We repeated the previous analysis at redshift z = 1. 
Fig.® displays these results in the same format as 
Fig.®. As before, for large values of the smoothing (e.g. 
Rs = 10 h~^ Mpc), the second-order correction is neg¬ 
ligible, while it becomes important for smaller smooth¬ 
ing scales when the local density around halos is high. 
In this case, choosing a smaller smoothing scale (e.g. 
Rs = 5 Mpc) and adding the second-order correction 
to Eq.® gives a better description of the displacement 
field (Eq.®). We note that at z = I the second-order 
correction is more important than at z = 0. This is 
due to the choice of halos we consider. At z = 1, ha¬ 
los with M > 10^^h~^MQ are more biased and belong 
to exponential tail of the halo mass function . Hence 
the local overdensity around those halos is larger, leading 
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FIG. 7: Same as Fig.®, for halos at redshift 2: = 1 in the 
mass range M > 


smoothing scales, the second-order correction becomes 
more important. 

Overall, at z = 1 the analytic approximations for the 
displacements (Eq.®) are in better agreement with the 
true displacements we measure in the simulations, by a 
factor which depends on the environment. For instance 
at the average density (E3, black cross), the standard 
deviation of the residual is ~ 2.bh~^ Mpc at z = 0 and 
^ 1.6/1“^ Mpc at z = 1. 


IV. EFFECT ON THE BARYON ACOUSTIC 
CORRELATION FUNCTION PEAK 

In this Section we study how the reconstruction of 
the baryon acoustic peak depends on the accuracy with 
which the displacement field can be determined. In par¬ 
ticular, we investigate if the “sharpness” of the recon¬ 
structed acoustic peak depends on local overdensity, and 
if this effect may be used to obtain improved accuracy in 
the resulting standard ruler measurement. 


A. The WizCOLA N-body simulations 



(3(Rs=5,7.5,10 Mpc.h-‘ ) 


FIG. 8: Same as Fig.®, for halos at redshift z = 1 in the 
mass range M > Mq. 


to a higher second-order correction of the displacement 
field. 

For halos with mass M > at z = 1, the 

error in the displacement is shown in Fig.®. For Rs = 
10 h~^ Mpc, the second-order correction is negligible and 
the error in the displacement is similar to that found 
for halos with M > 10^^h~^M q at z = 0. For smaller 


In order to obtain an accurate covariance matrix for 
our BAO fits, we require a much larger suite of N-body 
simulations than we used in the previous sections. We 
therefore analyzed the “WizCOLA” simulations ('[sijjl. 
which were produced in order to extend the Comoving 
Lagrangian Acceleration (COLA) technique [dj to lower- 
mass haloes to enable reconstruction of the baryon acous¬ 
tic peak in the WiggleZ Dark Energy Survey (Ei). The 
simulations are generated within a 600 h~^ Mpc box us¬ 
ing the “WMAP5” cosmological model. In particular, we 
used 1000 halo catalogues output at z = 0 in real-space, 
using halos with mass M > 10^^ h~^MQ as above, which 
yields a catalogue with number density ^ 5 x 10“^ 
Mpc“^, characteristic of current and future large-scale 
structure surveys at high redshift. 

We applied the density-field reconstruction method to 
these simulation boxes using standard methods described 
by, e.g. (Tsl - f^ . The process involves applying the dis¬ 
placements to both the data and a set of random points, 
which we generated uniformly within each simulation box 
with a number density 10 times greater than that of the 
data itself. The post-reconstruction correlation function 
is then measured using the displaced data and random 
points. When determining the displacement field we as¬ 
sumed a linear bias factor b ^ 1.4, fixed using the large- 
scale halo power spectrum. 


B. The correlation functions measured in different 
environments 

As in section UlI Cl we divided the simulation into dif¬ 
ferent density regimes by smoothing the density field 






traced by the halos. Fig.® displays the probability den¬ 
sity function n(i5, i?s) of the local overdensity 6 around 
halos (red histogram) and random tracers (black his¬ 
togram) at z = 0. The top panel is for a smoothing scale 
Rs = 7.5 h~^ Mpc, while the lower panel shows the PDF 
for Rs = 10h~^ Mpc. The vertical lines show the mean 
of the two distributions; for random tracers the mean is 
always zero, while for halos the mean of the distribution 
increases for smaller smoothing scales due to clustering. 

We split our sample into different environment bins 
using the local density around halos. For the purposes of 
this section we use 5 different environments. We defined 
the bin divisions by measuring the mean 5 and standard 
deviation as of the density values for each halo. In this 
case, environment corresponds to the average density 
and the central density values of the other environments 
are set using 

i5 = 5 -I- (* — S)as. (10) 

We also tested defining the bins based on the local den¬ 
sity around random tracers, which produced qualitatively 
similar results. 

We measured auto- and cross-correlation functions be¬ 
tween different environments using the Landy-Szalay es¬ 
timator (Hi). The cross-correlation function between 
environments {Ei,Ej) is therefore 

^ _ EEij TijiiTifij DRij Tisii DRji Tisij ^ ^ (11) 

RRij RRij RRij 

where DDij are the pair counts of halos between envi¬ 
ronments (i,j) noi is the total number of halos in the 
environment i, RRij are the pair counts of random trac¬ 
ers, DRij the cross-pair counts, and nai the total number 
of random tracers in environment i. 

We determined the covariance matrix of the cross¬ 
correlation function between two separation bins (r^, n) 
by averaging over the ensemble of N-body realizations: 

Ck,i = ^„ (0) - (6,(fc))(6j(0) (12) 

The error in the measurement of ^ij (jk) is then yjCk,k- 

In Fig. (HOD we plot the auto-correlation functions ^ii{r) 
measured in the 5 different environments after recon¬ 
struction. The black line (and dots) correspond to E 3 , 
the red lines (and dots) correspond to overdense envi¬ 
ronments E 4 , E^ while the blue lines (and dots) are for 
underdense environments Ei,E 2 . In order to produce 
a clearer visualization for the purposes of this plot, we 
stacked the simulation boxes in groups of 10, such that 
the error in the measurements is then determined from 
100 realizations of 10 stacked boxes. Fig. ra clearly 
displays a hierarchy where the baryon acoustic peak is 
sharper for the underdense regions and becomes more 
blurred for the high density regions. Interestingly, the 
most under-dense environment produces a “bell shape” 
around the acoustic peak. The cross-correlation func¬ 
tions display behaviour intermediate between the corre¬ 
sponding auto-correlation functions. 
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FIG. 9: The probability density function for the smoothed 
overdensity 5 around halos and random tracers measured at 
z = 0. The top panel shows this PDF for Rs = 7.5 h~^ Mpc, 
while the lower panel is for Rs = 10Mpc. 


In order to further understand how the correlation 
functions are affected by the reconstruction method, 
we plot in Fig. m the full correlation function before 
and after reconstruction. We see that reconstruction uni¬ 
formly improves the sharpness of the acoustic peak in all 
environments. 

The increased sharpness of the post-reconstruction 
baryon acoustic peak in underdense environments could 
potentially be used to improve the accuracy of the stan¬ 
dard ruler by upweighting those environments. However, 
the fluctuating weight could also serve to increase the co- 
variance in our measurements. We now explore the trade¬ 
off between these effects. We note that using density- 
dependent weights will always result in a biased estimator 
for the underlying correlation function. However, given 
that the amplitude and shape of the correlation function 
are marginalized over in standard methods of fitting the 
preferred scale, this bias is not necessarily problematic if 
the acoustic peak itself is not shifted. 


C. Total correlation function and new estimators 


The total correlation function ^ can be expressed as 
a function of the environmental correlation functions 
Ea. dTTl) . We begin with the Landy-Szalay estimator for 

c, 


£R<_2—^ + i 

RR njj RR 


(13) 
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FIG. 10: The auto-correlation functions measured in the dif¬ 
ferent environments, from bottom to top the measurements 
correspond to Cn,^55,^44,Css,C 22 . We use Rs = 10h~^ Mpc 
for the smoothing scale. 



r [Mpc.h^'] 


FIG. 11: Correlation functions measured in the different envi¬ 
ronments for Rs = 10 h~^ Mpc before (dashed lines) and after 
reconstruction (solid lines and dots). From bottom to top the 
measurements correspond to environments Ei, E 5 , E 4 , Es, E 2 . 


where DD, RR and DR are the total pair counts (halo- 
halo, random-random tracer, halo-random tracer). The 
total number of halos (nu), random tracers (n/j) and 
cross-product {njyn) appearing in this equation satisfy 

and = 

’’T'Di'riRj, where the sum is over all environments (and 
is written in this manner to allow us to generalize the re¬ 
lation as developed below). The conservation of the total 
number of pairs can be expressed as 

dd = J2dd.j ( 14 ) 

i,3 


Substituting Ea. dTTl) into Eo. lfOl) . we can express C as 
function of the Cy as 


f = !^J_ -2 

^ n^RR^ 




DR 

RR 


(15) 


where 


DDij = RR 




DR, 


nmnoj 

y Isy 4 ; 

nmriRj 


npi 
' nRi 


-DR 




riRj 


(16) 


1 . Illustrative example 

As a first illustration of the effect of assigning a dif¬ 
ferent weight to different environments, we re-write Cy 
in Eq. (IT5)l as Wij ^ij , where Wij = y/uHwJ is the weight 
assigned to each cross-correlation function in terms of 
weights Wi defined for each environment Ei. In the case 
where Wi = 1, we recover the original correlation func¬ 
tion. For the purposes of this example we consider a 
simpler split into 3 environments. 

In Fig. dT^ . we show the linear matter correlation func¬ 
tion rescaled by the linear bias (blue dashed line) and 
the non-linear matter power spectrum measured in the 
simulations (orange dashed line). The black line (and 
dots) shows the reconstructed correlation function using 
Rs = 10 h~^ Mpc with no weighting. For comparison 
we plot the weighted correlation function for two spe¬ 
cial cases (blue and red), where we choose wi = 2 , W 2 = 
2, W 3 = 0.5 for the blue line (upweighting low-density 
environments) and wi = 0,^2 = 0.5, rca = 2 for the 
red line (upweighting high-density environments). We 
can see that weighting the underdense environment (Ei) 
more than the overdense environment {E 3 ) allows us to 
recover the linear (sharpest) amplitude of the BAO peak. 


2. General linear weighting 


Hence we have obtained an expression for the total cor- We now construct a general estimator for the weighted 
relation function in terms of correlation function, insensitive to the absolute value of 
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FIG. 12: The linear correlation function (blue dashed line) 
and measured correlation function at z = 0 (orange dashed 
line). The standard reconstructed correlation function is 
shown by the black line, while the blue and red lines are par¬ 
ticular cases of our new estimate for the reconstructed corre¬ 
lation functions weighting by environment. 


FIG. 13: The linear correlation function (blue dashed line) 
and measured unweighted correlation function at z = 0 (or¬ 
ange dashed line). The solid lines show the weighted correla¬ 
tion function for x = —1, —0.5, 0, 0.5,1, which correspond to 
the blue, dark blue, black, orange and red lines, respectively. 
We do not show the errors in the measurement for clarity of 
the Figure. 


the weights, as 


where 


and 


- 


Sij i'^ij RRij + Pij) 
RRij 


nR 

riD 


n-Di n-Dj 

nRi URj 


(17) 


(18) 



( """I 

' (RDj 



\rdJ 

\nRj 


nm 


DRji + 

( riDi riDj 

riRi 

riDR^ 

\ nRi URj 




(19) 


and we re-defined njy = J 2 ij'^ij’^^Dinoj, n^R = 
J 2 ^ 3 '^ij'^R^'^Rj and = J^ijR’^jRDiRRj. This new 

estimator has the following desirable properties: ^ 

for Wij = constant, and = ^ki when Wij = SikSji. 


density environment, and a; is a variable that varies in 
the range x = [—1,1] such that for x = —1, the lowest- 
density environment is weighted with the highest ampli¬ 
tude = 2 and the highest environment is weighted 

with = 0. 

In Fig. dT^ we can see the resulting weighted correla¬ 
tion functions for various values of x, as well as the linear 
and non-linear correlation functions. We used 5 different 
environments and Rs = 10 h~^ Mpc. We observe that 
the hierarchy of the is not as simple as we show in 
Fig. dig. With this parametrization, the standard recon¬ 
struction (black line, x = 0) is closer to the linear corre¬ 
lation function than the measurement obtained assigning 
higher weight to the underdense environments (blue and 
light-blue lines) owing to the amplitude bias introduced 
by the weighting. However, the sharpness of the acoustic 
peak is nonetheless increased by such a weighting scheme. 
In the next section we consider the implications for the 
recovered standard ruler scale. 


To assign a value for each Wij = ^WiWj in Eg. (1171) we 
used the following parametrization: 


D. Baryon acoustic peak fitting 


Wi = 1 + X ( — - (20) 

\ ^max ^ / 

where i corresponds to the environment for which the lo¬ 
cal density 5{Rs) = <5, *max = 2i — fmin is the highest 


In order to determine the accuracy of the recovered dis¬ 
tance scale, we fit a simple baryon acoustic peak model to 
the galaxy correlation functions measured from each N- 
body realization, using the distribution of measurements 
across the realizations as the covariance matrix. Our 















fiducial correlation function model is based on a Fourier- 
transform of a model linear power spectrum P{k): 
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&d(s) 


1 

27r^ 


dkk^P{k) 


sin (fcs) 
ks 


( 21 ) 


We constructed the model power spectrum using a trans¬ 
fer function from the Eisenstein & Hu (1998) fitting for¬ 
mulae, assuming the fiducial cosmological parameters of 
the WiZCOLA simulations. We then fit for the scale dis¬ 
tortion parameter a, marginalizing over a normalization 
factor b^-. 


Cmod(s) = &^^fid(as) (22) 


This approximate model is adequate for our purposes 
of exploring the relative accuracy of determining a; we 
verified that a series of model extensions, such as using 
a model power spectrum from CAMB, or marginalizing 
over an additive scale-free polynomial or a free damping 
parameter, or using the mock mean correlation function 
as a template, did not qualitatively change our conclu¬ 
sions. 

We quantified the standard ruler performance by the 
standard deviation of the best-fitting values of a across 
the realizations, which we denote by a a. For this partic¬ 
ular test we used the 1000 single-box COLA realizations 
and 12 environments, checking that our conclusions were 
not sensitive to these choices. 

In order to test any improvement brought by the new 
concepts developed in the previous sections, we defined 
as a reference level the results obtained by a standard 
reconstruction method using the Zel’dovitch approxima¬ 
tion, smoothing length Rs = 10 h~^ Mpc, and no weight¬ 
ing (x = 0). Similar assumptions are used in standard 
implementations such as f [isl-l^i. Using the measure¬ 
ments of tTa, we computed the relative difference Acto, 
between the reference model and the results obtained (a) 
changing the smoothing length, (b) adding the 2LPT cor¬ 
rection and (c) introducing weighting of environments as 
a function of x. 

In Fig. dMl) we display our results for this analysis, in 
the form of the relative difference Actq, as function of 
the weighting parameter x defined by Eq. (1^ . The solid 
line uses the reconstruction method with the Zel’dovitch 
approximation and Rs = lQh~^ Mpc, the blue short 
dashed line shows the result for 2LPT and Rs = 7.5 
Mpc while the red dashed line shows 2LPT and Rs = 
5 / 1-1 

We find that the environmental-dependent weighting 
produces a small but measurable improvement in ctq, 
(~ 3%), with best performance usually produced in the 
range — 1 < x < —0.5. The implementation of 2LPT for 
computing the displacement field also produces a bene¬ 
fit as judged by CTo,, such that the total improvement is 
- 8 %. 



FIG. 14: The accuracy of the recovery of the standard ruler af¬ 
ter density-field reconstruction, relative to a reference imple¬ 
mentation, as a function of the parameter x in Eq. (Uni) which 
controls the relative weight assigned to underdense and over- 
dense regions. The reference case is the ZA with Rs = 10 h~^ 
Mpc and no weighting. 


V. DISCUSSION AND CONCLUSION 

We summarize the conclusions of our study as follows: 

• The displacement of proto-halos can be predicted 
with high accuracy by the ZA using their pecu¬ 
liar velocities. In this case, Eq. © evaluated at 
Rs = Rl provides a very good description of 
the displacement field with accuracy ~ 0 .bh~^ 
Mpc, independently of the halo masses. The choice 
of the smoothing scale is important, since there is 
a unique scale for which the approximation gives 
the best agreement with the exact final positions of 
the halos. The error in each component of iF fol¬ 
lows a Gaussian distribution centered on zero with 
a variance sensitive to the smoothing scale and the 
environment of the halos. 

• When reconstruction of the displacement field is 
performed using the z = 0 halo distribution in real 
space, we likewise find that the accuracy depends 
significantly on the smoothing scale and environ¬ 
ment. We extended the ZA to 2LPT and estab¬ 
lished a link between the error in predicting the 
displacement field and the local density around ha¬ 
los. 

• Based on this result, we showed that after applying 
reconstruction, the baryon acoustic peak is sharper 
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in the correlation function of low-density environ¬ 
ments. We hence built a new estimator of the cor¬ 
relation function, constructed by weighting a set 
of auto- and cross-correlation functions measured 
between different environments. 

• Fitting a BAO model to the results, we found a 
small but measurable improvement of ^ 8% in de¬ 
termining the standard ruler scale through a com¬ 
bination of using 2LPT to find the displacement 
field and weighting the environmental correlation 
functions. Further improvements may be possible 
using more sophisticated weighting schemes. 

The fact that the reconstruction of the displacement 
field using the halo density distribution at 2 ; = 0 is ~ 5 
times less accurate than the ZA applied to the initial ve¬ 
locities of proto-halos (independently of the orders con¬ 
sidered in the LPT approximation) implies that the mea¬ 
surement of the error contains non-linear 

information which could potentially be sensitive to dif¬ 
ferent cosmological models, especially in high-density en¬ 
vironments. The evolution in time of the environmental 
correlation functions should carry this non-linear infor¬ 


mation and can be use to test non-standard cosmologies. 

For instance in Fig. dm, the correlation function be¬ 
fore reconstruction in high-density environments has a 
smeared BAO peak compared to that of the lowest den¬ 
sity environment. The ratio of the BAO peak widths in 
those two environments should be directly proportional 
to the degree of late time non-linear interactions. For 
non-standard cosmologies, this ratio might be different 
(e.g., fewer non-linear interactions for early dark energy 
models). We will investigate these topics further in fu¬ 
ture work. 


VI. ACKNOWLEDGMENTS 

We warmly thank Y. Rasera for providing easy access 
to the DEUS simulations [3^. We are also grateful to 
Jun Koda for generating the WizCOLA simulation suite 
used in this investigation. Part of this research was con¬ 
ducted by the Australian Research Council Centre of Ex¬ 
cellence for All-sky Astrophysics (CAASTRO), through 
project number CE110001020. We also acknowledge sup¬ 
port from the DIM ACAV of the Region Ile-de-France. 


[1] T. Matsubara, Physical Review D, vol. 78, Issue 8, id. 
083519 

[2] Padmanabhan, Nikhil; White, Martin; Cohn, J. D., 
Physical Review D, vol. 79, Issue 6, id. 063523 

[3] Carlson, Jordan; Reid, Beth; White, Martin, Monthly 
Notices of the Royal Astronomical Society, Volume 429, 
Issue 2, p.1674-1685 

[4] K. C. Chan Physical Review D, Volume 89, Issue 8, 
id.083515 

[5] Eisenstein, Daniel J.; Seo, Hee-Jong; Sirko, Edwin; 
Spergel, David N., The Astrophysical Journal, Volume 
664, Issue 2, pp. 675-679. 

[6] Peebles, P. J. E.; Yu, J. T., Astrophysical Journal, vol. 
162, p.815 

[7] Eisenstein, Daniel J.; Hu, Wayne, The Astrophysical 
Journal, Volume 496, Issue 2, pp. 605-614. 

[8] Blake, Chris; Glazebrook, Karl, The Astrophysical Jour¬ 
nal, Volume 594, Issue 2, pp. 665-673 

[9] Seo, Hee-Jong; Eisenstein, Daniel J., The Astrophysical 
Journal, Volume 598, Issue 2, pp. 720-740 

[10] Seo, Hee-Jong; Eisenstein, Daniel J., The Astrophysical 
Journal, Volume 665, Issue 1, pp. 14-24. 

[11] Seo, Hee-Jong; Siegel, Ethan R.; Eisenstein, Daniel J.; 
White, Martin, The Astrophysical Journal, Volume 686, 
Issue 1, pp. 13-24 

[12] Smith R.E., Scoccimarro R., Sheth R., 2008, PhRvD, 77, 
43525 

[13] Angulo, R. E.; Baugh, C. M.; Frenk, C. S.; Lacey, C. 
G., Monthly Notices of the Royal Astronomical Society, 
Volume 383, Issue 2, pp. 755-776 

[14] Crocce, Martin; Scoccimarro, Roman, Physical Review 
D, vol. 77, Issue 2, id. 023533 

[15] Sanchez, Ariel G.; Cole, Shaun, Monthly Notices of the 


Royal Astronomical Society, Volume 385, Issue 2, pp. 
830-840 

[16] Y. B. Zeldovich, Astron. Astrophys. 5, 84 (1970). 

[17] White, Martin, Monthly Notices of the Royal Astronom¬ 
ical Society, Volume 439, Issue 4, p.3630-3640 

[18] Padmanabhan, Nikhil; Xu, Xiaoying; Eisenstein, Daniel 
J.; Scalzo, Richard; Cuesta, Antonio J.; Mehta, Kushal 
T.; Kazin, Eyal, Monthly Notices of the Royal Astronom¬ 
ical Society, Volume 427, Issue 3, pp. 2132-2145 

[19] Anderson, Lauren et ah. Monthly Notices of the Royal 
Astronomical Society, Volume 441, Issue 1, p.24-62 

[20] Kazin, Eyal A. et ah. Monthly Notices of the Royal As¬ 
tronomical Society, Volume 441, Issue 4, p.3524-3542 

[21] Noh, Yookyung; White, Martin; Padmanabhan, Nikhil, 
Physical Review D, vol. 80, Issue 12, id. 123501 

[22] McCullagh, Nuala; Szalay, Alexander S., The Astrophys¬ 
ical Journal, Volume 752, Issue 1, article id. 21, 6 pp. 
( 2012 ) 

[23] Tassev, Svetlin; Zaldarriaga, Matias, Journal of Cosmol¬ 
ogy and Astroparticle Physics, Issue 04, id. 013 (2012) 

[24] Burden, Angela; Percival, Will J.; Hewlett, Cullan, 
2015arXivl50402591B 

[25] Seo, Hee-Jong; Eckel, Jonathan; Eisenstein, Daniel J.; 
Mehta, Kushal; Metchnik, Marc; Padmanabhan, Nikhil; 
Pinto, Phillip; Takahashi, Ryuichi; White, Martin; Xu, 
Xiaoying, The Astrophysical Journal, Volume 720, Issue 
2, pp. 1650-1667 (2010) 

[26] Mehta, Kushal T.; Seo, Hee-Jong; Eckel, Jonathan; 
Eisenstein, Daniel J.; Metchnik, Marc; Pinto, Philip; Xu, 
Xiaoying, The Astrophysical Journal, Volume 734, Issue 
2, article id. 94, 12 pp. (2011) 

[27] Burden, A.; Percival, W. J.; Manera, M.; Cuesta, Anto¬ 
nio J.; Vargas Magana, Mariana; Ho, Shirley, Monthly 



13 


Notices of the Royal Astronomical Society, Volume 445, 
Issue 3, p.3152-3168 

[28] Vargas-Magana, Mariana et al.. Monthly Notices of the 
Royal Astronomical Society, Volume 445, Issue 1, p.2-28 

[29] Landy, Stephen D.; Szalay, Alexander S., Astrophysical 
Journal, Part 1 (ISSN 0004-637X), vol. 412, no. 1, p. 
64-71 

[30] White, Martin, Monthly Notices of the Royal Astronom¬ 
ical Society, Volume 450, Issue 4, p.3822-3828 

[31] Koda, J., et ah, 2015, in preparation 

[32] Alimi, J.-M.; Fzfa, A.; Boucher, V.; Rasera, Y.; Courtin, 
J.; Corasaniti, P.-S., Monthly Notices of the Royal As¬ 
tronomical Society, Volume 401, Issue 2, pp. 775-790. 

[33] http;//www.deus-consortium.org 

[34] Courtin, J.; Rasera, Y.; Alimi, J.-M.; Corasaniti, P.-S.; 
Boucher, V.; Fzfa, A., Monthly Notices of the Royal As¬ 
tronomical Society, Volume 410, Issue 3, pp. 1911-1931 

[35] Rasera, Y.; Alimi, J.-M.; Courtin, J.; Roy, F.; Corasan¬ 
iti, P.-S.; Fzfa, A.; Boucher, V., INVISIBLE UNIVERSE: 
Proceedings of the Conference. AIP Conference Proceed¬ 


ings, Volume 1241, pp. 1134-1139 (2010) 

[36] Teyssier, R.,Astronomy and Astrophysics, v.385, p.337- 
364 (2002) 

[37] Komatsu E. et ah,The Astrophysical Journal Supple¬ 
ment, Volume 180, Issue 2, pp. 330-376 (2009). 

[38] Spergel, D. N., et ah, 2007, ApJS, 170, 377 

[39] Rasera Y. et al., 2014, Monthly Notices of the Royal As¬ 
tronomical Society, Volume 440, Issue 2, p.1420-1434 

[40] Blot, L.; Corasaniti, P. S.; Alimi, J.-M.; Reverdy, V.; 
Rasera, Y., Monthly Notices of the Royal Astronomical 
Society, Volume 446, Issue 2, p.1756-1764 

[41] , Journal of Cosmology and Astroparticle Physics, Issue 
10, article id. 077, pp. (2014) 

[42] M.Kopp, C. Uhlemann and I. Achitouv, In prep 

[43] Achitouv, I.; Wagner, C.; Weller, J.; Rasera, Y., Journal 
of Cosmology and Astroparticle Physics, Issue 10, article 
id. 077, pp. (2014) 

[44] Tassev et al. 2013 JCAP 0636 


